arXiv: 1502.07013v2 [math.AP] 17 Sep 2015 


ENERGY SCALING LAW FOR THE REGULAR CONE 


HEINER OLBERMANN 


Abstract. We consider a thin elastic sheet in the shape of a disk whose reference 
metric is that of a singular cone. I.e., the reference metric is flat away from the center 
and has a defect there. We define a geometrically fully nonlinear free elastic energy, 
and investigate the scaling behavior of this energy as the thickness h tends to 0. We 
work with two simplifying assumptions: Firstly, we think of the deformed sheet as an 
immersed 2-dimensional Riemannian manifold in Euclidean 3-space and assume that the 
exponential map at the origin (the center of the sheet) supplies a coordinate chart for 
the whole manifold. Secondly, the energy functional penalizes the difference between the 
induced metric and the reference metric in L°° (instead of, as is usual, in L^). Under 
these assumptions, we show that the elastic energy per unit thickness of the regular cone 
in the leading order of h is given by C*h^ \ log h\, where the value of C* is given explicitly. 


1. Introduction 

1.1. Statement of the result. The present contribution is concerned with the shape 
that one obtains when removing a sector from a thin elastic sheet in the shape of a disc, 
and gluing the edges of the cut back together. We will call this setup the “regular cone”. 
It has previously been considered, e.g., to describe the buckling behavior of membrane net¬ 
works near disclinations, see 3^. In this reference, one can find ansatz-based discussions 
and numerical simulations of the regular cone. For more recent numerical simulations, 
see (^ . In ^|, minimizers of the free elastic energy have been described in some detail 
under the assumption of small deflections and rotational symmetry. The latter is a partic¬ 
ularly strong assumption on the shape of (almost-) minimizers. In the present work, this 
assumption will be removed. 

In a simplified setting, the elastic energy of a thin sheet consists of two terms, “mem¬ 
brane” and “bending” energy; the former is non-convex, and the latter can be viewed as 
a singular perturbation. By the non-convexity of the energy, there can be many local 
minimizers. One approach to a rigorous study of the problem, in the context of nonlinear 
elasticity first used in pfli . 21 1 , is to establish upper and lower bounds for the minimal 


energy in terms of the singular perturbation parameter. This is the approach we will take 
in the present paper. The perturbation parameter will be the sheet thickness h. 

We will define a free elastic energy for the regular cone and construct upper and lower 
bounds for the free elastic energy. The upper and lower bound will be identical to leading 
order in h. Our method of proof is based on viewing the elastic sheet as an immersed 
Riemannian manifold and analyzing intrinsically defined geometric objects, such as the 
Gauss curvature of the manifold. For our method to work, however, we are led to consider 
a non-standard free energy functional, and to restrict the set of allowed geometries of the 
deformed sheet. 
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The idea is to consider integrals of the Gauss curvature as the fundamental objects that are 
controlled by both the membrane and the bending energy (in different function spaces). 
The integral of Gauss curvature over some suitably smooth set U determines the (oriented) 
volume of the push forward of U under the Gauss map, i^{U). Via the isoperimetric in¬ 
equality on the sphere, this can be translated into a lower bound on on ||Z)z^||ii(g( 7 ), which 
will be enough to obtain lower bounds for the bending energy. 

For this idea to work, we need to be able to control the Gauss curvature by the membrane 
energy in some suitable function space. It is a well known fact from Riemannian geometry 
(already known to Gauss himself) that the Gauss curvature does indeed fully determine 
the metric (see Lemma[T]). However, this dependence only becomes explicit in a certain set 
of coordinates, the so called geodesic polar coordinates. In order to fully control the metric 
via the Gauss curvature on the whole manifold, we will have to assume that the geodesic 
polar coordinates supply a coordinate chart for the whole manifold. This is equivalent 
with requiring that there exist no no conjugate points to the origin (which is, roughly 
speaking, the same as requiring the non-existence of geodesics emanating from the origin 
that intersect elsewhere). 

To make this assumption precise, we introduce the notation Bp := {x G : \x\ < p} 
for p > 0, and we let y G C'°°(i?i;M^) be an immersion, i.e., rankDy = 2 everywhere in 
Bi. This map induces a metric on Hi, where y* denotes the pull-back under y, 

and denotes the Euclidean metric on M^. Now {Bi,y*e^^^) is a Riemannian manifold. 
We denote the tangent space at x G Hi by T^Bi. With this terminology, our working 
assumption can be formulated as follows: 

Assumption 1. For the Riemannian manifold (Hi,y*e^), the exponential map at the 
origin expQ supplies a diffeomorphism of some subset Ny of ToHi to Hi. 

Let A denote the set of immersions y G C'°°(Hi;M^) that satisfy Assumption [TJ 
We write gy = y*e^^^ for the induced metric, and set the reference metric go to be that of 
a singular cone, 

go = dp® dp + mlp'^d'd ® di? , (1) 

where p,^ are the polar coordinates on Hi and mo G (0,1). Furthermore, we denote by 
i^y the surface normal diy A d2y/\diy A d2y\. With this notation in place, we are ready to 
dehne the free elastic energy : C'°°(Hi;M^) —)• M by 

(y) =-^membrane + -^bending 

= \\9y - 9o\\Loa(^Bi\Bh) + ^^II-^^2/IIl2(Si) ' 

(For a precise definition of the norm on the space of metrics that we are using, see Section 
H) Furthermore we dehne the constant 

C* = 27r (l — mo) . 

In this contribution, we prove 
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Theorem 1. There exists a constant C > 0 depending only on tuq with the following 
property: For all h < 

C*\\ogh\ - ^log|log/i| - C <h-^ inf If°{y) < C*\logh\+C. 

2 yeA 

Of course, the choice of I'ff as the membrane energy requires further comment. 

Being able to use the geodesic polar coordinates globally by our Assumption (TJ it still 
remains the problem that the elastic energy should be given in the coordinates of the 
reference configuration. We will need to show that smallness of the membrane term 
implies that these two sets of coordinates differ only by little (in an L°° sense), in order 
to be able to translate energy estimates in one set of coordinates into the other. 
Unfortunately, we have not been able to do so using a more standard free energy like 

h{y) = WOy - 9 o\\l'2{Bi) + h‘^\\Dl^\\L'2{Bi) ■ ( 3 ) 

(For a justihcation of the latter, starting from 3-dimensional nonlinear elasticity, see 
e.g. (I^.l The reason why we have to work with instead of is that in the con¬ 
struction of the lower bound, we will have to assume that the metric error gy — g^ is small 
in and not only in Lf (cf. Remark [T]). 

Furthermore, in the definition of the membrane energy F^^^brane = Wdy “ 9o\\L^{Bi\By), 
we excluded a ball of radius h from the domain so that our upper bound construction (see 
Lemma [2]), that we believe to capture the qualitative shape of the regular cone, satisfies 
the “right” energy bound C*h?\ log/i| -|- Chf (see also Remark [2]). 


Next, we would like to comment on Assumption [TJ On the one hand, there exists a 
large set of deformations that satishes the assumption. On the other hand, note that the 
non-existence of conjugate points to the origin is a strong assumption in the following 
sense. Let K denote the Gauss curvature and dA the volume element of {Bi,gy). The 
integral |Ar|dA is the geometric quantity that is controlled by the bending energy. 
On a given geodesic, there will exist a conjugate point if enough (positive) curvature is 
concentrated on that geodesic. This can be achieved by a modification of the deformation 
that raises |iF|dA by an arbitrarily small amount. Hence, the existence of conjugate 
points is not penalized by the energy (at least not sufficiently to assume that low energy 
maps satisfy Assumption [1] in general). 

While so far we have only been able to carry out the interpolation technique for the Gauss 
curvature in geodesic polar coordinates, we believe that the feasibility of this approach 
should in principle not depend on the choice of coordinates. The reason for this hope is 
that the objects in our interpolation inequalities (i.e., the metric and the Gauss curvature) 
have a coordinate independent geometric meaning. 


1.2. Scientific context. Starting in the 1990’s, there has been a lot of activity in the 
physics and engineering community devoted to the description of energy focusing in thin 
elastic sheets. The regions where the elastic energy is focused are approximately straight 
folds, approximately conical singularities, or a (possibly very complex) network consisting 
of both. The first instance of a scaling law for the elastic energy in such a settin g, in the 
sense of an asymptotic limit for vanishing sheet thickness h, can be found in |27l . 1^ . 
The first detailed analysis of the energy scaling for a “minimal ridge” in the Foppl- 
von Karman setting is due to Lobkovsky [2^. Single ridges were further investigated 

37l |. The vertices where elastic energy focuses were studied in a 
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simplified setting (the so-called d-cone) in JIMJ . There also exist earlier works on the 
buckling transitions for spherical shells [^, sjsjthat display energy focusing in certain 
regimes. A thorough overview over this area and many more references can be found in 
the very recommendable review article by Witten [^ . 

The mathematical works closest to the present one in spirit and content are [^, [l3, • 

In [I^, it has been shown that the elastic energy per unit thickness of a “single fold” 
scales with , building on results from 371. In [§, 321, the following has been proved: 
Consider an elastic sheet in the shape of a disc, and a conical deformation with the apex 
of the cone at the center of the sheet. Now restricting to non-singular configurations of 
the sheet that agree with the singular one on the boundary and at the center, the elastic 
energy per unit thickness scales with h‘^\ log/i|. 

Mathematically, the papers [^, 12, treat variational problems for energy functionals 
that have a membrane and a bending term, with certain Dirichlet boundary conditions. 
These boundary conditions are chosen such that there exists a unique (Lipschitz) isometric 
immersion that satisfies them. Additionally, this map is singular in the sense that it has 
infinite bending energy. (In the sequel, we will call such boundary conditions “tensile”.) 
Deviations from the isometric immersion that satisfies the tensile boundary conditions are 
penalized by the membrane energy. The optimal balance between membrane and bending 
energy yields the lower bound in the elastic energy. 

Obviously, this method of proof relies heavily on the choice of boundary conditions. The 
method completely brakes down in the current setting. 


To the best of our knowledge, the use of tensile boundary conditions has so far been 
the only tool to prove lower bounds in variational problems for thin elastic sheets that dis¬ 
play energy focusing in ridges or vertices. We speculate that the interpolation techniques 
for the Gauss curvature that we use here can be a tool for proving such lower bounds when 
tensile boundary conditions are not present. This would be particularly interesting in the 
so-called crumpling problem. The latter consists in determining the minimal amount of 
energy necessary to ht an elastic sheet into a container whose diameter is smaller than the 
diameter of the sheet. 

Of course, to make the speculation about wider applicability of our techniques more sub¬ 
stantive, we would have to get rid of the technical assumptions made in the present work. 
This will be the topic of future research. 


1.3. The relation to the (7^ “ isometric immersion problem. In this subsection, we 
want to comment on the relation between the present results and the isometric immersion 
problem with low regul a,rity . 

By the results by Nash |33| and Kuiper (^ . every short immersion i/ : t/ —)■ of {U^go) 
can be approximated arbitrarily well in the C^-norm by C^-isometric immersions. This is 
an instance of Gromov’s famous h-principle [T^. Looking at the elastic energy functional 
/“, the Nash-Kuiper construction yields maps with arbitrarily small membrane energy 
G^-close to a given short immersion y. Hence, there is an abundance of maps with arbi¬ 
trarily small membrane energy. On a heuristic level, the fundamental difficulty that we 
encounter here, and that is shared by the present problem with the crumpling problem, is 
how to show that all of these many degrees of freedom have large bending energy. 

In sharp contrast to the setting, there is the classical rigidity in the Weyl problem. 
This result states that any isometric immersion of a given Riemannian manifold {S‘^,g) in 
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is rigid, i.e., unique up to a rigid motion. 

The striking contrast between the results for C^- and isometric immersions natural^ 
leads to the question how the situation looks like for (7^’“ isometric immersions. In [7|, 
Borisov announced that if g is analytic, the h-principle holds for isometric immersions in 
for a < A proof can be found in [g]. Another result by Borisov [^, 
states that for a > |, any isometric embedding y G of a C^-metric g on 

with positive Gauss curvature is rigid. 

For our purposes, the latter is relevant. The proof of the statement consists in show¬ 
ing that the embedded manifold M is of bounded extrinsic curvature, i.e., 

N 

sup {Y,^\^yiEi)) : N€N, 

i=l 

a collection of closed disjoint subsets of m| < oo , 

where denotes the two-dimensional Hausdorff measure. This requirement is sufficient 
to obtain the rigidity result, using classical results by Pogorelov [s^. We see that quanti¬ 
tative control over extrinsic curvature (the left hand side above) eliminates the possibility 
of constructions in the style of Nash and Kuiper. Indeed, extrinsic curvature diverges in 
the Nash-Kuiper approximation scheme. Since we will be dealing with smooth maps here, 
the extrinsic curvature is just the L^-norm of the Gauss curvature, \Ky\dAy, which is 
controlled by the bending energy. Heuristically, this supports the view that gaining con¬ 
trol over the Gauss curvature is also the natural way of dealing with variational problems 
that allow for short maps. 

1.4. Plan of the paper. This paper is structured as follows: In Section [SJ we dehne 
our notation and make some preliminary observations and definitions. In Section [3l we 
establish L°°-estimates for the exponential map, that in the sequel will allow us to translate 
bounds on the elastic energy into estimates of the metric and the Gauss curvature in the 
geodesic polar coordinates, see Lemma [3l In Section [H we use these results to interpolate 
between the metric and the Gauss curvature, as in the ansatz that we have proposed 
in this introduction. One more ingredient is needed to carry out the ansatz, namely an 
appropriate isoperimetric inequality on the sphere, which is proved in Section [5j In Section 
E we combine these items to prove our main result. 

Acknowledgments. The author would like to thank Stefan Muller for helpful discus¬ 
sions. 


2. Notation and preliminaries 

For a manifold M and x G M, T^M denotes the tangent space at x, and T*M its dual. 
For a hnite dimensional vector space V, the space of alternating fc-linear forms M 

is denoted by A^V. A /c-form on a manifold M is a map that associates to every x G M 
an element of A^T^M. If M is a subset of M^, the gradient Df of a function f : M ^ R™ 
is dehned via duality to df with the standard Euclidean metric, 

Df={d,f,d2ff, 
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even if M is equipped with a Riemannian metric that is different from the Euclidean one. 
The /c-dimensional Hausdorff and Lebesgue measures will be denoted by 1-L^, respec¬ 
tively. The oscillation of a function u : U —)• M is denoted by o^c u = sup[/ u — infj/ u. The 

space of Lipschitz functions from U C M"’ to some metric space V is denoted by Lip ([/;!/); 
if R = M, we simply write Lip (U). 

Unless stated otherwise, the symbol “C” will be used as follows: The statement “/ < 
is short-hand for “There exists a constant C > 0, that only depends on mg, with the 
property / < Cg”. The value of C may change from line to line. Sometimes it will be 
convenient to be able to refer to the same constant later on in the text; in this case, we 
denote it by Cj, i = 1, 2,..., and it will be fixed. 

We now want to make precise the definition of the polar coordinates p, "d, the reference 
metric go, and of the norm on the space of metrics, that have already been used in the 
introduction. 

On Bi, we introduce polar coordinates p, 'd by 

:Ri^(0,l) 

X l-^(|x|, x/|x|). 

The symbols dp, will at the same time denote partial differentiation with respect to p, d 
as well as the vector fields Bi —)• TRi, defined by their actions on functions / : —)• M, 

{dpf) (x) = 1^ • Df{x), {di)f) {x) = • Df{x), 

where x-^ = (x 2 , —xi)^ for x = (xi,X 2 )'^. 

Now we fix 0 < mg < 1 and define the reference metric gfg : —)• T*Bi (g) T*Bi by ([2]). 

Next we need to specify the norm on the space of metrics used in the definition of the free 
elastic energy Q. For x G \{0}, we define a norm on inner products p : T^Bi x T^Bi —)• 
M on the tangent space at x: For 

p =piAp (g) d/3 -|- P2P{(^P ® dd -|- di? (g) dp) + pgp^dd (g) dd 

we set 

IIpII^ = Pi + 2p2 + pi • 

This makes || • || well defined since dp{x), pdd(x) span T*Bi for every x G \ {0}. Then 
we define the L°°-norm of a metric g on A <Z Bi hy 

IIs'IIl-(A) =esssup,^g^||5r(x)|| . 

The objects Ny, (expg)y ,Vy as well as the maps that we will introduce in the following 
all depend on y ~ for notational convenience we will not indicate this dependence anymore 
from now on. 

Next we introduce radial coordinates {r,p) with respect to C TgRi — M^, 

(r, ip) ■. Bi\ {0} -^(0, oo) X 

X (|x|, x/ |x|) , 

where x = (expg)“^(x). 

The vector fields dr,dy, are defined on N analogously to the definition of dp,d^ on Bi 
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above; and using the push forward expg, we will view them as vector fields on Bi from 
now on. 

By Assumption [U the map {p{x),'&{x)) i-A {r{x),(p{x)) is regular on its domain and hence 
we may write 

=r( 

where T : Bi\ {0} —)> is defined by 

dpr p~^d^r \ 
rdpip rp~^d^ip ) 

Note that 





rdr A d(^ = det Tpdp A di? . 

By Assumption [H rdr A d(p is a nowhere vanishing two-form on Bi \ {0}, and hence 


det r > 0 . (5) 

We also introduce the following notation for the inverse of B: Let V \ B\ \ {0} —)• be 

defined by 

^ ^ ( drp r-^d^p \ 

Y por'd pr~^d^'d j 

An obvious consequence of this definitions is 

r(x) = for all x G i?i \ {0} . 

Note that in the relation above, L and have the same argument, since we understand 
all partial derivatives as functions on the manifold Bi. Hence, 

_ 1 ( pr-^d^'d -pOr^ \ 

detf V OrP J ' 


Let iL ; Hi —)• M denote the Gauss curvature, and dA the volume 2-form of (Hi, i^). 

We are now going to define three functions Q,Cl,G : Hi —>■ M, and we will make these 
definitions with respect to the geodesic polar coordinates {r,ip). By Assumption [H this 
will make them well dehned as functions on Hi. 

We set 

H(0, (/?) = H(0, (f) = 1 — G(0, (/?) = 0 for all ip £ 

and 


H(r, p) 

n{r,p) 


rr 

■■ / K{s,p)dA^^^^){dr, 
Jo 

1 r 

/ H(s,(/?)ds 

^ Jo 


G=i-n. 


dy,)ds 


( 7 ) 
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3. Passage to geodesic polar coordinates 

The proof of Theorem [1] will be an application of the ansatz that we proposed in the in¬ 
troduction - however, to make the interpolation argument work, we need to pass from the 
polar coordinates in the reference configuration (p, i?) to the geodesic polar coordinates 
(r, (/?). This is what we will do in the present section. 


The starting point is to note that in the (r, (/?) coordinates, the metric takes a partic¬ 
ularly simple form. The following lemma can be found as Proposition 3 and Remark 1 in 
Chapter 4 of (I^ . or in Chapter 3 of For the convenience of the reader, we include 
the proof here. In this proof, we denote the Levi-Civita connection of g by V. 


Lemma 1. The metric g : Bi ^ T*Bi ^ T*Bi is given by 

g = dr iSi dr + G'^r'^dtp (g) dip . 


Proof. First, we have g{dr,dr) = 1 by the definition of the exponential map. Secondly, 
g and its first derivatives are identical to the Euclidean metric at the origin, by basic 
properties of the exponential map. Hence 


lha g\ (d^^dr) = 0 for all ip. 


Next, we compute the derivative of g{d^,dr) along geodesics, 


drg{d^, dr) =g{Vrd^, dr) + g{d^,Vrdr 


=0 


=g{V^dr,dr) 

=~dipg{dr, dr) 
=0 


( 8 ) 


In the second equality above, we used that the Lie brackets of the coordinate vector fields 
vanish: 

[dr, - Vy,dr = 0 , 

and that dr is a geodesic vector field, Vrdr = 0. The initial conditions together with (jSj) 
prove that g{d^,dr) = 0. It remains to show that = G^. 

For vector fields X,Y let R{X,Y) = [Vx,Vy] — X[x,y] denote the Riemann curvature 
tensor. We compute 

Xld^=Vr(y^dr) + Vr[d^,dr] 

=R{dr, d^)dr + V^Vrdr 

= -Kd^. 

The value of dA evaluated on the vector fields dr, d^p and the first derivatives of these 
values w.r.t. r at the origin are the same as in Euclidean space, by the defining property 
of the exponential map. Hence 

limd^|(.^((9,.,5^) =0 

r—>-0 ' ' 

lima^dAb ) ((9^,(9^) =1 

r—>-0 ^ 
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We compute the second derivative of d^(9r.,9(^) along the radial curves parametrized by 
r: 


drdA{dr,d^) =dA{dr,Vld^)) 

= - KdA{dr,d^) 


( 10 ) 


Together with the initial conditions ([9]), this defines an initial value problem, and thus 
dA{dr,d^) satisfies 



By the definition of the volume form, 



□ 


The following lemma serves two purposes: On the one hand, it supplies the upper bound 
in the statement of Theorem [1] On the other hand, it assures that the assumptions in the 
subsequent lemmas make sense. 

Lemma 2. We have 

inf (y) <C*h‘^\logh\+ Ch^ . 
y&A 

Proof. We define the following il-dependent orthonormal frame in 

Cp = (cos-djSini?, 0 ), = (—sini?, cosi9,0), 62 = ( 0 , 0 , 1 ) 

Further, we set Cmo := iR-oCp + y^l — mge 2 . Let if G C'°°(M+) such that tpip) = 0 for 
p < 1/2, il){p) = 1 for yo > 1 and \ilj'{p)\ < 4, |'0"(/3)| < 8 for all r. We claim that the 
upper bound is satisfied by the map defined by (using polar coordinates on Bi) 


y{p,'&) := 'ilj{p/h)pemo + (1 “ '^{p/h))pep . 


We calculate 



+ (1 - (1 - 1710 ) 1 ^) e,9 (g) e,9. 

Hence the pullback of the Euclidean metric in under y is given by 



(-0 + + (1 - V’ 




dp (g) dp 


+ (1 — (1 — p^di} (g) di? . 


(11) 


The surface normal is given by 


{^ + W emo + (l-V’-FV’0 

fip) 
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where 


fip) = y (l - ml) (ij + + (mo + 1 - V' - ■ 

We compute the derivative of the normal, 


DUy = 




(gmo - 


fip) 

ih'^' + fimo + (1 - - W) g 

fip)^ 


(2 - mo) (^-0 + - (1 + m-o) 


( 12 ) 


1 (V’ + fV’O v/i - mi 




■e^ 0$ 6,9 


P fip) 

From m and m, and using the properties of 'i/’j we see that 

IlS'y - fi'olP = 0 forp>/i, 

< C for p < h 




= h?p ^(1 — mg) for p > h 


Thus we get 


I'hiv)= sup \\gy-gof + h?l \Dvy\^ + h^ 

B\\Bh J Bh 

p2tt p1 

' / di? / 

Jo J h 


\DvA 


Bh 
(1 -m^) 


’Bi\Bh 




□ 


Bi\Bh 

<Ch^ + 

=C*h^\\ogh\+Ch^. 

This proves the present lemma. 

In the following Lemma, we use the smallness of the membrane term \\g — to 

prove a certain “rigidity” of the exponential map expg. Namely, we prove that circles in 
N are approximately mapped to circles in Bi in an L°° sense. This allows us to pass from 
the radial coordinate p in the reference configuration to the radial coordinate r in A^. 

Lemma 3. Let y ^ A with I^iy) < 2C*h‘^\ log/i|. Setting tq := supg^^^ have 

sup |r(p,-!?) — ro — p| < Ci/il log. (13) 

2h<p<l 

i?eSi 

Proof. Let dg,dgg denote the distance functions with respect to g,go respectively, i.e., 
dgix,x') = g\-yit)il'it)A'it))^^‘^dt : 7 G Lip ([0,1]; Bi), 7(0) = x, 7(1) = x'j . 

The minimum is achieved by a geodesic connecting x and x'. The function dg^ is defined 
analogously. 

Let X £ Bi \ Bh- Recall that r{x) is the distance from the origin with respect to the 
induced metric g, 

r{x) = dg{x,t)). 

We claim that 

osc r < Ch , 
dB2h 
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Indeed, let x,x' £ dB 2 h- There exists a curve 7 £ Lip ([0,1]; \ Bh) with 7(0) = x, 

7(1) = x', and goil'(t),'j' (t)) < Ch'^ for all t £ [0,1]. Hence 

Jo Jo 

<Ch. 

Hence, by the triangle inequality for dg, \dg{x,d) — dg(x',0)| < dg{x,x') < Ch, which 
proves (H- 

Now let X £ Hi \ B 2 h, and 7 a geodesic connecting 0 and x with 7(0) = 0 and 7(1) = x. 
There exists some to G (0,1) such that l7(to)l = 2/i and |7(t)| > 2h for all t £ (to, !]• We 
write xo = 7(to), and note ||x| — (igp(x,xo)| < Ch. 

Now we have 


|r(x) - ro - |x|| <|dg(x,xo) - dga{x,xo) + dg{0,xo) - ro| + Ch 
<\dg{x,xo) -dga{x,xo)\ +Ch, 

where we have used (|14p in the second inequality. 

Let 7o be a curve connecting x and xq with 7o(0) = xo, 7o(l) = x and 

dg^{x,xo)= [ 5ol7oW(7c)(0>7o(0)^''^dt. 


(15) 


Then 


dg{x,xo)= [ 5l7oW(7oW>7o(^))^'^^dt 
Jo 

<Jo ((1 + <^/i|log^l^^^)5ol7oh)(7oW>7o 

<{1 + Ch\ logti|^/^)(tgo(x,xo). 


1/2 


dt 


(16) 


On the other hand, for any curve 7 connecting x and xo with 7o(0) = xq, 7o(1) = x, we 
have 


^ fi'l77)(7'W,7'W)^'^^dt (^(1-C'h|log/i|^/2)5rol7(7(7'W,7' 

>(1 + Ch\ logh|^/^)dgo(x,xo). 


1/2 


dt 


Hence we also get dg{x,xo) > (1 — Ch\ log/i|^/^)ctgp(x, xq). Using this last inequality and 
()16p in ()15p . the claim of the lemma is proved. □ 

Remark 1. If one works with the energy functional Ih coming from the Kirchhoff-Love 
ansatz instead, it is possible to show an analogous claim if one additionally assumes that 

11^90^ - ep||L2 < C'llfi'-5'o||L2 , (17) 

where Vg^ denotes the covariant derivative w.r.t. g^, and Cp denotes the unit vector in p 
direction in the reference configuration. 

We believe that dm) does hold true, since 


\\g - gof < Udist^ Grvl if ) ’ (pointwise) , 
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cf. equation (EU) in the proof of Lemma below. This looks quite similar to the situation 
in the Geometric Rigidity Theorem by Friesecke, James and Muller 0 /. However, we 
have not been able to adapt these ideas to deduce (HZD from m- 

If one is able to show (dzi), then one can deduce more or less the same statement as 
in Lemma 0 through a combination of arguments as in the proof of Morrey’s inequality 
11 /11(^0,i-«/p < C'll/lliyi.p and AssumptionUl Then one would also be able to show the same 
statement for Ih as in Theorem U\ (with different coefficients in front of the lower order 
terms). 


4. Interpolation between metric and Gauss curvature 


Now we are in position to follow the ansatz from the introduction and interpolate between 
the metric and the Gauss curvature in the (r, (p) coordinates. The interpolation will be of 
the standard type 


(see e.g. 


171, Theorem 7.28), where 


(19) 


fir) = —[ dipG{r,ip) 
mo J51 


linear terms in r . 


Since G is defined via the double integral of the Gauss curvature (cf. jT])), the second 
derivative dff{r) can be bounded in L^ by the bending energy. 


In order to establish smallness of ||/||li from the smallness of the metric error g — go, 
we will need the following lemma. 


Lemma 4. There exists a constant C = C{mo) with the following property: 
(i) For every y ^ A, 


/ Gdetr 
V mo 



< C\\g-go 


(pointwise), 


and hence for every y G A with Iff{y) < 2C*h?‘\ log/i|, 


Gdetr 

mo 


< Ch\\ogh\^/‘^ . 


(ii) For every y G A with If^{y) < 2C*h‘^\ log/ij, every R with 2h < R < 1, and every 
dGS\ 



\dpr{p, 'll) — l\dp < Ch\ log . 


Proof. First we introduce the notation 


mo 


1 

0 



G 


1 0 
0 G 


By the definition of go and Lemma [U 

50 = (dp, pdi?) (g) Tho ^ ^ , g = (dr, rd(p) (g) 

By definition of the matrix norm || • || (cf. section [2|), we get 


dr \ 
rd(p J 
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where on the right hand side, | • | denotes the usual Euclidean matrix norm. By multiplying 
— ?71 q from the left and right with we alter the norm of this expression at 

most by a factor tuq, and get 

\Q'^Q-ldf <C\\g-go f , 

where 

Q = GVfriQ^ , 

and the constant on the right hand side only depends on mo. Now we write down the 
spectral decomposition of the symmetric matrix Q^Q, 

Q^Q = DR, 

where D = diag((ii,d2) is some diagonal matrix and R G SO{2). By ||A|| = ||ii^74i?|| for 
all A G and R G SO{2), and detQ^Q = did 2 we easily deduce 

iVdetQ^Q-lp < |Q^Q-Id|2 

and thus 

. ( 20 ) 

mo ) 

This proves the hrst claim of (i), and the second claim of (i) follows trivially. 

Next we note 

|Q^Q-Id|2 =|(x/^-Id)(x/Q^ + Id)|2 

>\4WQ - Id|^ = dist2(Q,0(2)) 


where the inequality holds by \JQ^Q + Id > Id in the sense of positive dehnite matrices. 
Hence we have shown 

dist2(Q,0(2)) <0||5-<7of . 

In fact, since det Q > 0, we even have 

dist2(Q,50(2)) <0||<7-5of . (21) 


Explicitly, Q reads 


/ dpr ttiq \ 

Y rGdpp rGmQ^p~^d^Lp ) 


( 22 ) 


We introduce the notation (•)+ = max(-,0), (•)_ = max(—,0). From (I2TI1 and (l22]) . we 
get in particular 

{dpT- 1 )_^ < G\\g- goW < Gh\\ogh\^G _ 

As in the proof of Lemma O we set ro := supao r and note that oscr < Gh. Obviously, 

B2h 

we have dp{r — p — tq) = dpV — 1. Hence, for every R with 2h < R < 1, and fixed i? G 
(the arguments p, being omitted in the notation), 



rR 

l)_dp < sup |r — p — ro| + osc r + / {dpV — 

9B2h J2h 


l)+dp 


<Ch\logh\^/^ 

13 
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rR 


Hence we get 

[ \dpr-l\dp= f {{dpr-l)^ + {dpr-l)_)dp 
J2h J2h 

<Ch\\ogh\^/‘^. 

This proves (ii) and completes the proof of the present lemma. □ 

We will need to distinguish balls in Bi and in N C TqBi, and hence we write 

Br ={x G Bi : p{x) < ii} 

Br ={x G Bi : r{x) < ii} . 

Also, for a measurable set V C M, let 

/C(V) = KdA. 

In the statement of the following proposition let, as in the proof of Lemma 01 

ro = sup r{2h, ?i). 

As a further piece of notation, we set 

r* := 1 — 2h + vq — Cih\ log , 

where Ci is the constant from the right hand side of (|13p . By this choice, Br* C Bi. 

Proposition 1. Let y £ A with Ih{y) < 2C*h?‘\ log h\. Then for R G [Cih\ log h\^/'^, 1 — 
2h - Cih\\ogh\^/% 

r2h+R 

/ dp |/C(i?p) — 27r(l — mo)| < log. 

J2h 

Proof. Let 

p = p{r,(p), T} = '&{r,ip) 

be understood as functions of the coordinates (r, <p). By Lemma El 


sup |p(r, (p) + ro — r| < C/il log. 

rG[ro,r*], 


(23) 


Now define 


/ : 


G 


r Ir d<p- I — 27r(r — tq) 

JdBr . 


We may rewrite 


/(r) = r [ 

Jd 


dBr \^o r 


dp A 


Note that 


f'(r) = ( — f dpfl ) — 27r = (— f iLdA ) — 27r 

KmoJaBr J \moJBr J 

f''{r) =— [ dtpKGr . 

1^0 JdB^ 


14 
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By the Gauss equation, 


K = 


det Du^Dy 


det Dy'^Dy ’ 

and since for all x € Bi, the image of Di^lx is contained in the image of Dy\x, we have as 
a consequence 

\K\dA = det Du'^Dy ^ Vdet Du^Dudx . 

x/det DyTDy 

Now we can estimate \\f"\\L^(ro,r*) as follows, 


ll/"llLi(.o.r-*) = r drdy^\KGr\ 
J rn 

\K\dA 


L 

L 


Br* \BrQ 


dxV det Du^Di 


(24) 


B^* \BrQ 

< [ dxlDi^l"^ 

Jbi 

<C\ logh \, 

where we have used the upper bound for the bending energy in the last inequality. 

For the application of the standard interpolation inequality (I19D . we need to estimate the 
norm of /. More precisely, we will estimate ||/||Li(ro,ro+R)- 

rro+R rro+R 




dip 


rro+R 

< / rdr 

Jro 

rro+R 

< / rdr 

Jro 

rro+R 

< / rdrdy? det F 

Jro 


\1 


r 

r — tq r 


det TdpV 


+ 


r-To 


\mQ r p 

V mo p 

' GdetF ^ 
mo 

r - ro - p 


dpv 


(25) 


\dpr - 1| + 


P 


where we used ([6]) in the second estimate to obtain the relation d^iJ = detTrp~^dpr. We 
will estimate the three terms in the integrand on the right hand side separately. Before 
we do so, note that by Lemma O 


Bro+R B {Q,R + 2h + Ch\ log 


C\BR+ro) <c(R + h\\ogh\^/‘^y 


<CR^. 

15 
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Now we estimate the terms on the right hand side in (j25p . Firstly, 


pro+R 

GdetF 


GdetF 

/ r dr dp det i' 


< 


Iro 

mo 


mo 




C-^{Bro+R) 


<Ch\\ogh\^/'^R^ , 

where in the second inequality, we have used the assumption I^{y) < 2C*h‘^ \ log/i| and 
Lemmad] (i). Secondly, setting Rh = R + 2h + Cih\ log 


\dpr - 1 

•Rh 


pro+R 

/ rdrdp det F 

o 

1 

/ ro 

p 


< sup 

r>ro 


r — Vq 


p 

rRh 


f pdp f di? \dpr — 1| 
J2h isi 


/ 

J2h 


r — 1| 


<CRh 
<CRh\ log/i|^/2^ 

where in the last inequality, we have used Lemma d] (ii). Thirdly, 

r-ro+R 


Jtq 


rdrdip det F 


1 

o 

1 

Cl. 

VI 

di? 

1 

o 

1 

p 

J2h J, 

Ji 

p 


rRh 

=C / dp\r-ro- p\ 

J2h 

<CRh sup |r - ro - p\ 
<Ch\logh\^/^R. 

Inserting the preceding estimates into (f2^ , we get 

Using (IMl) and (f27)l in the interpolation inequality (fT9]l we obtain 

\\IC{Bp) - 27r(l - mo)|Ul(.o,ro+R) = ll/'llLl(r-o,r-o+R) < ChV2i?l/2| log/,|3/4 


(26) 


(27) 

(28) 


Comparing this last estimate with the statement of the present proposition, we see that 
what remains to be done is a change of domain for the integration of KdA, from Bf = 
{x £ Bi : r{x) < f} to Bp = {x £ Bi : p{x) < p}. We have 


r2h+R r2h+R / 

/ d/9|/C(Sp) - 27r(l - mo)| < / dp JC{Bp) - JC{Bp+ro- 2 h) 

J2h J2h \ 


(29) 


+ 


p+ro-2h) 27r(l — mo) 


In the following estimate we will use the notation AAB = {A\ B) U {B \ A) for the 
symmetric difference of sets A, B, and the characteristic function of a set A is denoted by 

16 
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XA- The first term on the right hand side of (I29p can be estimated as follows, 


r2h 

J2h 


2h+R 


dp 


— }C{Bp^ro-2h) 


p2h+R r 

/ - 

J2h J E 


< 


I2h 

r2h+R / r 

J2h \J 

<Ch\logh\^/^. 


^p-\-rQ — 2h p 


\K\Grdrdip 


) \K\Grdrdp 


(30) 

where in the first estimate, we have used the definition of /C, in the second line we have 
used Fubini’s Theorem, and in the last inequality, we have used Lemma [3] and (|24p . The 
second term on the right hand side of (I29|) can be estimated by a trivial change of variables 
and 


r2h 

J2h 


2h+R 


rro+R 


dr 


JC { Br ) — 27r(l — mo) 


(31) 


nr 

dp K,{Bp+ro-2h) - 27r(l - mo) = / 

J ro 

Inserting ([3n|l and m into (1291) concludes the proof of the present proposition. □ 


5. An isoperimetric inequality on the sphere 


In the proof of our main theorem, we are going to need a certain isoperimetric inequality 
for the image of the Gauss map. We will need several items from the literature to prove 
our claim in Section [TH Lemma [7] below; they will be collected in Sections 15.1115.31 

5.1. The isoperimetric inequality for sets on the sphere. Let 5^ = {x G : I x\ = 

1}. We define F : [0,47r] —)• M by 

F{x) = yFKx — a^. 

Note that F is concave. The isoperimetric inequality for sets on the sphere is as follows: 
Theorem 2 ([3]). Let A d he open. Then 

n^{dA) >F{n^{A)). 


5.2. The Brouwer degree. The statement of the isoperimetric inequality will involve 
the Brouwer degree. We mention some basic facts about this object - for a more thorough 
exposition with proofs of the claims made here, see e.g. (I3 |. 

Let M be a paracompact oriented manifold of dimension n, and U a bounded subset of M”’. 
Further, let u G C°°{U; M). Assume that y £ M \ u{dU), and let p he a C°° n-form on 
M with support in the same connected component of M\u{dU) as y, such that p = 1. 
Then the degree is defined by 


deg(y, U, u) 



where u* is the pull-back by u. It can be shown that this definition is independent of the 
choice of p. Further, deg{-, U,u) is constant on connected components of M \ u{dU) and 
integer valued. Moreover, it is invariant under homotopies, i.e., given H G C°°([0, 1] x 
U, M) such that y ^ H{[0, l],dU), we have 


deg{y,U,H{0,-))=deg{y,U,H{l,-)). 

17 
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Using these facts, one can go on to define the degree for u £ C^{U, M) by approximation. 
Since the reader may be more familiar with a different way of defining the degree, we 
mention that for points y £ M \ u{dU) that satisfy deiDu{x) ^ 0 for all x £ rt“^({y}), it 
may be shown that 

deg(y, U, u) = E sgndet Du{x). 

xeu-^{{y}) 

This can also be used as a starting point to define the degree. 

If u £ C^{U\ M), and C^{dU) = 0, then it follows from the defining formula for the degree 
and approximation by smooth functions that 


!■ 


deg{y,U,u)iJ.= / u*(/r) 


(32) 


'U 


for any n-form y that can be written as 

fj-iy) = A • • • A d?/n 

in local coordinates {yi ,..., yn), with ip £ L°°{M). Assuming that M C with m > n, 
the above formula implies 

f I deg(y, [/, u)|d^”(y) < f \ del Du{x)\dC^{x). 

Jm Ju 

We conclude that deg(-, U,u) £ LJ{M) for u £ C^{U;M). If y in is an exact form, 
i.e., 

fi = duj 

for some n — 1 form cj on M, then 

u*(da;) = d {u*u) . 

If U has smooth boundary, this implies, by Stokes’ Theorem, 



U, u)dui = f d {u*ui) 

Ju 

L 


= U UJ . 

IdU 


(33) 


5.3. The space BV{S‘^). In the following, we are going to consider the space L^{S‘^) = 
L^{S ‘^where T-L^ is the 2-dimensional Hausdorff measure. For / £ L^(5^), let 

\\Df\\{S^):=mi!^liminil^jDMdn^: f, G Up (S^), f, ^ f G L\S^)^ . (34) 

Then BV{S'^), the space of functions of bounded variation on the 2-sphere, is the set of 
/ £ L\S^) with WDfWiS^) < oo (see P, |3^, where functions of bounded variation are 
defined for much more general measure spaces). 

We will need some basic facts about BV{S‘^), that we state in Lemmas MM below. We 
first need to collect some notation for the statement of the first part of Lemma [3 
For X G S‘^, the space TxS'^ is naturally equipped with the norm coming from the standard 
metric g on S^, 

= V9\x{v,v) for V G . 

18 
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The space /\^TxS'^ is equipped with the dual norm 

II«IIait ^52 = sup{|a(u)| : v G T^S'^, < 1} for a G ^^TxS‘^ . 

Now define C^(5^; A^TS^), the space of C^-one forms on 5^, to be the set of functions 
oj that associate to every x G an element of with the property that for any 

coordinate chart : 5^ D 17 —)• we have 

X i-A ('0*cu)|3;(ej) G C^{U) for i = 1, 2, 

where e* denotes the unit vector in i-th direction in On A^T5^), we introduce 

the norm || • ||o,i by setting 

lli^llo,! := sup ||a;|3;||/\i7’^52 . 

xGS^ 

Lemma 5. (i) For all f G 

WDfWiS^) = sup /(iu : a; G C\S^-, A^TS^), ||a;||o,i < 1 

(ii) If A C S'^ is open and has Lipschitz boundary, then 

\\DxA\\iS^)=nHdA), 

where XA denotes the characteristic function of the set A. 

Proof For all f G BV{R'^), 

||i7/|| (M^) = sup I f / div h dx :/i G (7^(1^^; M^), ||h||o < 1 

Ur2 

where || • || denotes the C^-norm. For a proof of this statement, see e.g. [^. The statement 
(i) follows using a smooth atlas on 5^ and a subordinate partition of unity. Statement (ii) 
follows from (i) and an application of the Gauss-Green Theorem. □ 


The following coarea formula for BV functions is taken from [30|], where it is proved for 
much more general measure spaces. 

Lemma 6 For all f G 

/ OO 

ds\\Dxif>s}\\{S^), 

-OO 

where X{f>s} denotes the characteristic function of the set {x G : /(x) > s}. 

Next, we claim that for u G C^{Br; S'^), the Brouwer degree deg{-, Br,u) is in BV{S‘^). 
Indeed, we first note that 'H'^{u{dBfi)) = 0 and hence deg{-, Br, u) is defined 7^^-almost 
everywhere on S‘^ and from (j32p . one can easily derive that it is in L^(S‘^). Using 


/ deg(-, Bp, u)da; = / u*da; 
752 JBp 

= [ d(M*a;) 


L 

I 

JdBo 


U U} , 
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where, in the second equality, we have use the fact that pull-back and exterior derivative 
commute, and in the last equality, we have used Stokes’ Theorem. Hence, by Lemma [5] 
(i), we have 

\\D{deg{-,Bp,u))\\iS‘^)= [ \Du\dn^ for all u € C\Bp; S^). (35) 

JdBp 


5.4. An isoperimetric inequality for Z- valued functions on S'^. As a last piece 
of notation before we make the main statement of the present section, we introduce the 
metric d47rZ on the real numbers given by 


diirzix) = dist(a;, dvrZ). 


Lemma 7. Let Br C v G C^{Bpt,] 5^) and 


Then 


/C = / deg{x, BR,u)d3i‘^{x). 
Js^ 



> F{d4^z{}C)) . 


(36) 


Proof. We set 

As := {x G 5^ : deg(x, Br, z^) > s} , 

and note that Ag is open and has Lipschitz boundary. Combining Lemmas [5] (ii) and [6] 
with equation (f35]l . we get 


By Theorem [2j we get 


Now 


[ \Dn\dn^ 

JdBn 


'H^{dAs)ds 


' —oo 
oo 


j=-oo 



oo 

> ^ F{n^{Aj)) 

j=-oo 


(37) 


}C=f deg{x,BR,i^)d'H‘^{x) 

Js^ 

= Y,'F\A,)-Y,'FHS^\Aj), 

j>0 j<0 


(38) 


where the right hand side is an absolutely convergent sum by deg{-, Br,^) G L^(5^). By 
the triangle inequality for d47rz. 


^47rZ (^) sE {Aj)) -I- d^T^zj—TL^\ Aj)) 

j>0 j<0 

= ^ d4Trz{'H^ {Aj)) . 

j€Z 
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By the concavity of F, we obtain 

F o d^TrziK.) < E F o d4^z{n\Aj)). 


Since G [0, Att] for all j and F = F o d^T^z on [0,47r], this can be rewritten as 

Fid^MlC)) <Y.F{n\Aj)). 

j£Z 

Combining this last inequality with (|37p proves the claim of the lemma. 

6. Proof of the main Theorem 


□ 


Proof of Theorem [H The upper bound has been proved in Lemma [2j Let y satisfy the 
upper bound. For any R G [Cih\ log 1 — 2h — Cih\ log we have 


r2h-\-R 

/ dp |/ C ( Bp )- 27r ( l - mo )| < Ch^/^R^/^\logh\^/^ 
J2h 


by Proposition [TJ We set F = F o dist47rz with F as in Lemma [71 Note that on 
Imo ■= ~ ^)) P is Lipschitz with a Lipschitz constant that only depends on 

772-0. Moreover, F is bounded by dvr. Hence, 

fcpTno)!® - 27r(l - ?72o)| if x G 1^0 
Ibvr^ else. 


(f{x)) -( T ( 27r ( l - mo )))' 
This implies 


< 


r2h 

J2h 


2h+R 


dp 


/ ^ \ 2 /_ \2 r2h-\-R 

( F (/ C ( Hp))j - ( T ( 27r(l - mo ))) <C J dp |/ C ( Bp ) - 27r(l - mo )| 

<C/2^/2^^/2nQg^|3/4_ 


(39) 


Now we may estimate the bending term, first using Jensen’s inequality: 

On the right hand side, we may apply Lemma [71 and obtain 

[ /■‘dp 

JBi 


-F{)C{Bp))^ . 

2vr Jo p 

We set hi = 2Cih\ log/i^/^, and choose J G N such that 

2'^hi < l-Ci/2|log/2|i/2 <2'^+l/li. 

Note that this choice implies in particular 

3 

log 2'^ > I log h\ — - log I log h\ — C. 
21 
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Using this in (Hfl . we get 



1 

> — 
“27r 

>C* 


hi 


P 


nnBp)f 


1 

log 2-' - — / ^ F{JC{Bp)f - F(27r(l - mo))' 

27r Jhi P 


>C* { I log ^1 “ 2 I ~ F ) — error term . 


(41) 


Letting Rj = 2^hi, the error term on the right hand side can be estimated using (I39p . 



P 


F{)C{Bp)y^ 


F(27r(l - mo))^ 



P 


nnBp)f 


F(27r(l - mo))^ 


< ^ R-^ChB'^{2Rj)F ^I log h\ 

j=0 

<C, 


where the last estimate is just the summation of a geometric series. This proves the 
theorem. □ 


Remark 2. It is apparent from the proof that we would have been able to prove the same 
lower bound (with different coefficients in front of the lower order terms) if instead of 
working with the membrane term 115—5011^00(5^^5^); we had worked with II5—5011500(5^^^5,); 

where h = Ch\ log hf for some p > 0. 
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